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The Bulge Radial Velocity Assay (BRAVA): I. Sample Selection 

and a Rotation Curve 
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ABSTRACT 



Results from the ongoing Bulge Radial Velocity Assay (BRAVA) are pre- 
sented. BRAVA uses M red giant stars, selected from the 2MASS catalog to lie 

CN \ within a bound of reddening corrected color and luminosity, as targets for the 

Cerro Tololo Inter-American Observatory 4-m Hydra multi-object spectrograph. 
Three years of observations investigate the kinematics of the Galactic bulge ma- 
jor (-10°</<+10°, h=-A°) and minor (-6°<6<+5°, -0.4°</<-0°) axes with 
~3300 radial velocities from 32 bulge fields and one disk field. We construct 
j^ \ a longitude- velocity plot for the bulge stars and find that, contrary to previous 

studies, the bulge does not rotate as a solid body; from — 4°</<+4° the rotation 
curve has a slope of roughly 100 km s~^ kpc^^ and flattens considerably at greater 

f->. ' /, reaching a maximum rotation of 75 km s~^. We compare our rotation curve 



and velocity dispersion profile to both the self-consistent model of Zhao (1996) 



C^ ' and to N-body models; neither fits both our observed rotation curve and velocity 

C^ ' dispersion profile. We place the bulge on the plot of (Vmax/o") vs. e and find 

^yTj ■ that the bulge lies near the oblate rotator line, and very close to the parameters 

O ■ of NGC 4565, an edge-on spiral galaxy with a bulge similar to that of the Milky 

Way. 
/\ • We find that our summed velocity distribution of bulge stars appears to be 

J-H ' _1 

c^ ■ sampled from a Gaussian distribution, with a = 116 ± 2 km s~ for our summed 

bulge fields. Furthermore, the high precision of our radial velocities (~5 km s~^) 
allows us to investigate hints of cold kinematic features that were seen in a 
number of the line-of-sight velocity distributions from our earlier observations. 



^Department of Physics and Astronomy, University of California, Los Angeles, CA 90095-1562 

^Visiting Astronomers, Cerro Tololo Inter-American Observatory. CTIO is operated by AURA, Inc. under 
contract to the National Science Foundation. 

■^ Cerro Tololo Inter- American Observatory, Casilla 603, La Serena, Chile 

^SUPA, School of Physics and Astronomy, University of St. Andrews, KY16 9SS, UK 



-2- 



In the past, optical radial velocity studies in the bulge have not emphasized high 
precision, because of the large velocity dispersion and the expectation that the 
short orbital periods would erase any cold structures in well under a gigayear. 
Thus our precision is sufficient to enable a search for more cold streams analogous 
to those associated with the Sagittarius dwarf spheroidal galaxy; some candidate 
cold features were seen and follow-up observations are reported herein. 

Subject headings: Galaxy: bulge — Galaxy: kinematics and dynamics — stars: 
kinematics — stars: late-type — techniques: radial velocities 



1. Introduction 

The status of the Milky Way bulge as a distinct stellar population had historically been 
well established by Baade's surveys of RR Lyrae stars, and the well known concentration 
of red giants toward the Sagittarius region (Blanco 1965). Furthermore, the image of the 
Galaxy produced by the COBE satellite (Weiland et al.l984; Dwek et al. 1995; Arendt et al. 
1998) depicted an unambiguous, asymmetric, peanut-shaped bulge. Blitz & Spergel (1991) 
modeled the 2.4/im balloon borne data of Matsumoto et al. (1982) as a bar which shows the 
hallmark thickening of the bulge at positive Galactic longitude that is seen in later studies 
(see e.g. Launhardt et al. 2002). While the discovery of a bar in the stellar distribution was 
somewhat of a surprise, gas motions in the inner 2 kpc have long been modeled as due to 
the potential of tilted bar (e.g. Liszt & Burton 1980). Optical surveys find roughly 1/3 of 
spirals are barred (Sellwood & Wilkinson 1993); this fraction rises to 60% or more in the 
infrared (Menendez-Delmestre et al. 2007), so it is perhaps not surprising that one should 
find evidence of bar structure in the Milky Way. In fact, the evidence for a triaxial (barred) 
bulge is convincing. The COBE 2/im light distribution is modeled as a bar oriented toward 
positive Galactic longitude (Dwek et al. 1995, Binney et al. 1997). Star counts of red 
clump stars (Stanek et al 1997; Babusiaux & Gilmore 2005) show a barlike structure, while 
microlensing events further show a central bar that is pointed roughly in the direction of the 
Sun (Figured]). 

There is also the issue of bar survival; Sheth et al. (2007) find that the bar fraction 
declines dramatically with redshift to half the Hubble time, with only the most massive 
galaxies hosting bars by 2; ~ 0.8. This fact is to be contrasted with growing evidence that 
the bulge is old from a chemical (Rich 1990; McWilliam & Rich 1994; Matteucci et al. 1999; 
Fulbright et al. 2007; Lecureur et al. 2007) and stellar populations (Ortolani et al. 1995; 
Kuijken & Rich 2002; Zoccali et al. 2003) perspective. There is also continuing debate over 
whether the Milky Way bulge/bar system is best categorized as a pseudobulge (Kormendy & 
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Kennicutt 2004). As evidence of a bar in the Milky Way has grown (see e.g. Gerhard et al. 
2002) it has become clear that the Galactic bulge offers a potential laboratory to investigate 
the dynamics of the nearest bar in detail with radial velocities and eventually with proper 
motions as well as abundances. This recent evidence has prompted renewed interest in the 
dynamical modeling of the Galactic bar/bulge system in order to determine, in part, the 
dynamical formation and evolution of galaxies e.g., the unsolved issue of whether a Cold 
Dark Matter halo cusp is compatible with the microlensing and rotation of the Milky Way 
bar (Zhao, Spergel, Rich 1995, Klypin, Zhao, Somerville 2002). 

However, Zhao's (1996) dynamical model paper notes the lack of comprehensive surveys 
as a problem in comparing the model against data; until recently there simply has not 
been enough data to properly constrain the few dynamical models of the bulge available. 
Indeed, many efforts have modeled the 2 micron light in projection (cf. Gerhard 2002), but 
surveys of the dynamics have generally either emphasized very rare evolved stars (e.g. SiO 
masers; Izumiura et al. 1995; OH/IR or IRAS stars; Menzies 1995) or use K/M giants. 
Previous studies are restricted to small isolated fields due to the problems of crowding, 
large and variable extinction, lack of wide field surveys in the red, and the presence of a 
contaminating disk population. Minniti (1992) used K giants, but at the time, uniform 
selection of candidates was an issue, and concerns were raised about disk contamination, 
especially for the field at (/,6) = (12°,3°). Consequently, past surveys have had insufficient 
statistics to properly probe the dynamics of fields at Galactic latitude | / 1 < 5° (see Table 3 
of Zhao 1996 for detailed summary) in order to provide a constraint for dynamical models. 

Partly in response to this problem, Beaulieau et al. (2000) undertook the first dynamical 
survey of the bulge aimed at isolating a uniform sample of kinematic probes while at the 
same time covering a wide area. Their study combines PNe from a new survey and those 
found in existing catalogs to span a wide range in Galactic longitude and to study both 
the Northern and Southern portions of the bulge. They argue that the easily observed PNe 
have no metallicity bias. However, PNe are short-lived and are therefore rare, and their 
population membership (and especially distances) are uncertain. Rich et al. (2007) noted 
(and we confirm) that the bulge PNe have lower velocity dispersion than other populations, 
consistent with a partial population of non-bulge members. 

The dynamical model for the bulge/bar has a number of important implications. Large 
samples of uniform radial velocity data are still of great value in constraining the bar versus 
axisymmetric models, and the nature of the orbit families supporting the bar. Further, the 
interpretation of the microlensing events in the bulge depends on the use of an accurate 
dynamical model (Han & Gould 2003). The recent discovery of planetary transit host stars 
in the bulge (Sahu et al. 2006) gives an additional incentive to improve our knowledge of 
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the bulge/bar model, as it is debated whether stars on non-circular orbits are incapable of 
bearing planets, if their orbits stray too near the Galactic center. The relatioinship to the 
disk, thick disk, and halo populations is of great interest. Only with a realistic dynamical 
model can the issue of long term stability be explored. There is also the issue of whether 
the bar formed due to the buckling of the disk (see e.g. Raha et al. 1991) or is a yet more 
ancient population. 

The availability of bulge photometry and high precision astrometry from the Two- 
Micron All-sky Survey (2MASS; Skrutskie et al. 2006) along with the Schlegel et al. (1998) 
all sky reddening map offered, for the first time, the possibility of easily selecting a large 
sample of kinematic probes from the very population responsible for the 2/xm radiation 
comprising the projected bulge/bar (Dwek et al. 1995; Launhardt et al. 2003; Figs. 1,2). 
The Bulge Radial Velocity Assay or BRAVA (Rich et al. 2007) was conceived as a survey 
of the line of sight velocity distribution of red giants across the bulge, to be compared with 
self-consistent dynamical models like that of Zhao (1996). We emphasize that while we 
report the first and second moments of the BRAVA fields in this paper, the ultimate aim is 
to study the Line of Sight Velocity Distribution (LOSVD) and to compare it with dynamical 
models. A further additional aim of BRAVA is to explore correlations between kinematics 
and abundances. Note that there is already a hint from one proper motion field study (Soto, 
Rich, & Kuijken 2007) that predominantly metal rich stars are observed to show a vertex 
deviation in their velocity ellipsoid and may therefore be predominantly responsible for the 
support of the bar. If [a/Fe] can be derived in addition to [Fe/H] , it may be possible to use 
the a— enhancement as a secondary criterion to distinguish bulge from thick disk stars, via 
chemistry (see e.g. Fulbright et al. 2007). 

Here we report results from a survey based on red giants, which comprise the bulk of 
the 2.4/im light of the bulge. M giants offer a dynamical probe that is populous, luminous, 
traces the light, and is easily utilized for velocity measurements. M giants are also extremely 
well studied (Frogel & Whitford 1987) and have been identified, and their giant branch 
characterized, over the whole of the bulge (Frogel et al. 1990), from —3° to —12° (420 to 
1675 pc, adopting Rq = S kpc). The issue of the abundance "bias" raised by Beaulieau et al. 
is essentially settled based on high dispersion infrared spectra (see Rich & Origlia 2005, and 
Rich, Origlia, & Valenti 2007). M giants follow roughly the abundance distribution of the 
K giants; given their significant contribution to the integrated light, we are confident that 
these stars are excellent dynamical probes. The M giants, being luminous, are also usable 
as probes even in regions of high extinction, and for exploring the fringes of the bulge, as 
Frogel et al. (1990) as well as the 2MASS catalog, illustrate. 

Early results from our survey are given in Rich et al. (2007a,b). Here, we report on the 
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first three- year results from BRAVA and also describe the sample selection and measurement 
uncertainties in detail. At this time, we present ~3300 spectra for which we have obtained a 
radial velocity precision of ~5 km s~^ for our most recent data. Future papers in this series 
will address the kinematics of the bulge at —8° and a detailed comparison with the Zhao 
(1996) model (Howard et al. 2008 in prep.). The long term plan is to make the entire survey, 
including the spectra, available on a public website. 

In the following section we discuss our sample selection. Section 3 discusses our ob- 
servations, while Section 4 discusses our data reduction and velocity calibration. Section 5 
considers the fields and the issue of reddening. Section 6 presents our results, while Section 
7 gives our conclusions. 



2. Sample Selection 

The first use of M giants to undertake a rapid survey of the velocity dispersion in 
Baade's Window (BW) was that of Mould (1983) who used the du Pont 2.5 m telescope 
at Las Campanas, and the photon-counting Shectograph, to obtain spectra of 50 stars and 
measure a velocity dispersion of a = 113 ± 11 km sec~^. The approach used in our present 
work was pioneered by Sharpies, Walker, & Cropper (1990, hereafter SWC), who used the 
AAT and a multifiber spectrograph to obtain spectra of a sample of 239 late-type M giants 
in BW. SWC found these M giants to be a kinematically distinct group in the bulge, and 
consequently excellent probes for studying the dynamics of the Milky Way bulge. Beaulieau 
et al. (2000) use the planetary nebulae (PNe) population of the bulge as dynamical tracers, 
using a sample of nearly 400 PNe in the extended region —30° < / < 30° and 3.3° < |&| < 
15° to investigate the kinematics. However, the brief lifetimes of PNe, uncertain population 
membership, and distance uncertainties make them a potentially problematic population 
to probe bulge kinematics. As a result, dynamical models of the bulge/bar system are, to 
date, relatively weakly constrained by radial velocity data. A better solution is to conduct a 
survey using red giant stars, following the SWC strategy, which are the most common type 
of luminous evolved star. We choose to work brighter than the red clump for two reasons. 
First, brighter giants can be studied in obscured fields, and in highly obscured fields, they 
can be studied in the infrared. Second, clump stars lying closer than the bulge tend to be 
both bluer and brighter, due to the proximity as well as the lower reddening. The RGB 
slants redward, the red edge of the RGB is comprised of the most metal rich red giants, 
largely uncontaminated by any other population. 

The most luminous red giants are cool enough to form titanium-oxide (TiO) molecules 
in their atmospheres, which is a signature of their spectral classification M. Because M giants 
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are extremely cool, their spectra are dominated by these TiO absorption lines, particularly 
towards the red end of the spectrum at 7055A, although other strong bandheads lie redwards 
of 8OOOA. Additionally, other spectral features such as the prominent calcium triplet at 
8498A, 8542A, and 8662A make M giants good candidates for radial velocity probes and allow 
for subsequent metallicity measurements. M giants, classified from low dispersion spectra, are 
known to be ubiquitous throughout the bulge (Blanco, McCarthy, & Blanco 1984; Frogel et 
al. 1990), and although faint at optical wavelengths due to the heavy extinction towards the 
Galactic center, TiO bands, and cool temperatures, they are relatively bright in the / band 
(~6500A to ~10000A). ^M^^^ gives us adequate astrometry to use the M giant population 
as dynamical probes of the Galactic bulge; these stars' relatively short lifetimes makes their 
numbers ~ 10 times less numerous than stars on the red giant branch and the clump, and 
consequently source confusion becomes a non-issue. M giants are also more luminous; even 
in the ~ 8OOOA / band where TiO bands are present, the M giant population is ~ 5 mag 
brighter than the red clump, another factor that suppresses crowding and complications 
compared to the red clump giants. This can easily be seen by comparing a 2MASS H-band 
image of Baade's Window with any optical image of that region. 

When choosing stars as dynamical probes, one must be careful to minimize foreground 
contamination by main sequence disk stars and giants in the near disk, as well as to minimize 
any metallicity selection bias. SWC observed Baade's Window and found that M giants with 
I < 11.8 have a lower velocity dispersion and thus were likely to be disk members, whereas 
the stars with I > 11.8 are likely bulge members, based on their kinematics. When the 
survey fields are examined in 2MASS, the Kvs. J- if color magnitude diagram (CMD) shows 
a clearly defined red giant branch (Figure [3]). Due to differential reddening towards the 
Galactic Bulge, the magnitude limit of / < 11.8 corresponds roughly to K < 8.2. Therefore, 
we adopt a range of 9.25> K >8.2 as the criterion for selecting 2MASS M giants; in principle, 
one can go even fainter, but the abundance of sources in the magnitude selection range is more 
than adequate, allows for less integration time, and avoids the red clump, an evolutionary 
path available only to stars with [Fe/H]> —1. 

Because the reddening varies greatly in the galactic bulge (1.5 < A„ < 5), a parallelogram- 
shaped region within the CMD is adjusted by eye in both magnitude and color to center 
on the expected locus of the distance and average reddening of each particular field, for 
which we used the extinction law calculator on the NED website (see Schlegel et al. 1998 for 
details). In fields with large extinction, the width of the parallelogram is widened to account 
for the large differential reddening. This initial selection method assures bulge membership 
and minimizes any metallicity bias. Our survey, to date, dynamically samples the Galactic 
bulge by observing red giants in fields at 1° intervals along the major axis at h =—4° along 
with the minor axis at /~0° (Figure [2]). 
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2.1. Disk Contamination 

The major issue in tlie selection of bulge samples is the purity of a bulge sample. The 
sightline crosses through the foreground disk, and these populations (mostly < 1 Gyr). Early 
samples were optically selected, and in those cases, there may be considerable uncertainty in 
assigning population membership. Sample contamination can arise from a number of sources. 
Thick disk stars with similar abundances, kinematics, age, and distance, are indistinguishable 
from the bulge. Such a population must be included in models of the total luminosity and 
dynamics. Thin disk populations are separable; they are younger and therefore, bluer. Most 
such contamination arises blueward of the most metal poor isochrone. But in highly reddened 
fields, differential reddening can scatter some thin disk stars into the selection zone of the 
bulge; dereddening has removed most of these. Foreground stars may overlay the selection 
region, by being closer and less reddened than the sample; this is a concern for the red clump, 
which is populous. We are 3 mag brighter than the red clump, so this is not an issue. The 
bar does have spatial depth, but all bar members will ultimately be modeled; kinematics as 
a function of photometric distance is beyond the scope of this paper. 

We have empirically controlled for a possible contaminating population by determining 
how purely Gaussian the velocity distributions are, and we have also compared the kinematics 
of subdivisions of the population by bright and faint, and blue and red (Section 6). The 
/ — V plot discussed in section 6 also shows no hints of extra subpopulations that might 
appear as linear features (a cold disk population). We have also studied a disk field at 
{l,b)= (—30°, —4°), finding its velocity dispersion a factor of 2 lower than any in our sample 
(Section 5). We are satisfied that our selection method yields a sample that is dominated by 
bulge/bar members and minimizes any disk contamination. Future extension of a study like 
this to include measurements of [Fe/H] and composition might be able to bring additional 
tools to bear in the assignment of population membership. 



3. Observations 

Observations began in 2005 and continued through 2007 (see Table [T]), yielding 32 
individual bulge fields and one disk field at {l,b) (—30°,— 4°). We use the Hydra multifiber 
bench spectrograph at the Cassegrain focus of the Cerro Tololo Inter- American Observatory 
(CTIO) Blanco 4-m telescope. On average, 103 stars are observed per field for a total 
of ~3400 stars observed. To take advantage of the red colors of M giants, we employ 
the KPGLD grating, blazed at 8500A giving a dispersion 0.45A per pixel with 2 pixel on 
chip binning, yielding an effective dispersion of 0.88A per pixel and a full spectral range of 
~1800A. In 2005, our central wavelength was ~7600A with an effective resolution of R~2800. 



However, in 2006 we adjusted redward with a central wavelength of ~7800A in order to 
observe the calcium triplet missed in 2005. In order to retain the TiO band at ~7050A, 
we were only able to observe the first two lines of the triplet. Despite this, very strong 
cross correlation peaks were obtained from the addition of the first two lines of the calcium 
triplet observed in 2006. In 2007 we adjusted further redward, with a central wavelength of 
~7900A. Additionally, in 2006/2007 we utilized the 200/im slit plate giving us an increase in 
resolution (R~4200). The light lost from the slit plate was determined to have little effect 
on our S/N due to the brightness of the sources. Typically, we obtained 3 exposures of each 
field at 600 seconds each, although a few fields required 3x900 second exposures, or longer, 
due to poor observing conditions. Each field has, on average, 105 successfully exposed M 
giants and ~20 dedicated fibers that are used to obtain a sky background spectrum. Table 
[2] shows a summary of our observations including the number of reliable velocities measured 
in each field. A sample selection of our stars from each year is shown in Figure HI along with 
a blue star to show regions of telluric absorption. The 2005 observations did not include any 
velocity standards, although we obtained 4 in 2006, and 2 in 2007 (Table [3]). 



4. Data Reduction and Velocity Calibration 

The raw frames are first trimmed and overscan corrected using the IRAF routine 'ccd- 
proc', and are bias corrected. Flat fielding, sky subtraction, throughput correction, scattered 
light subtraction, and wavelength calibration are all accomplished through the IRAF task 
'dohydra'. The spectra are binned to ~34.5 km s~^ pixel"^, co-added, and continuum nor- 
malized by a 2nd or 3rd order polynomial. Radial velocities are then measured using the 
IRAF cross correlation routine 'fxcor' utilizing a fourier filtering set to exclude both features 
exceeding 50 pixels (overall continuum) or smaller than 3 pixels (the spectral resolution). 
Regions of the spectrum contaminated by telluric features (such as the atmospheric A band 
at ~7600A) were omitted from the fitting, leaving roughly 60% of our total spectrum us- 
able for cross correlation; velocities are then determined by measuring a Gaussian fit to the 
strongest correlation peak (see fig [5]), and corrected to the heliocentric rest frame. 

All 4 velocity standards observed in 2006 are used in the 2006 cross-correlations. For 
the 22 2006 fields observed, 3 of the standards (HD 207076, HD 218541, and HD 203638) 
return individual stellar velocities that agree to better than 2 km s~^, on average, whereas 
the typical errors in stellar velocities reported by 'fxcor', as determined by the Tonry-Davis 
R- value, are of order of 5-10 km s~^. HD 177017 returns velocities that are offset from the 
aforementioned 3 standards by ~7 km s^^ on average with a standard deviation of 1 km s^^. 
Because velocity offsets found with HD 177017 were constant from field to field in the 22 
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fields observed in 2006, final 2006 velocities were calculated by applying a zero-point shift 
of the velocities returned by HD 177017, and then taking the weighted average of all four 
standards. 

Due to the different wavelength coverage in the 2005 and 2006 data sets, the 2005 
spectra have a smaller spectral range (~1450A) usable for cross correlation. Additionally, 
the lower resolution of the 2005 data combined with the loss of the calcium triplet results 
in lower quality cross correlation fits for the 2005 data. Errors in velocities as reported 
by 'fxcor' for 2005 are, on average 12-16 km s~^. The same velocity shifts applied to the 

2006 data are applied to the 2005 data, for consistency, and again, a final stellar velocity 
was computed from the weighted average. One field was observed in 2005 and 2006 at 
{l,b) = {Q°,—4°), with the same fiber configuration (i.e. the same stars in that field observed 
each year). To determine the consistency of our velocity results, a star by star comparison 
of stellar velocities of that field was conducted (Figure [6]). As can be seen, despite our lower 
resolution and less effective wavelength coverage in 2005, our stellar velocities show an offset 
of -2.7 km s~^, with an RMS scatter of 4.8 km s~^ about that mean. Since the offset is less 
than the RMS scatter, we consider the 2005/2006 data sets to be in good agreement, and 
adopt 5 km s~^ as our errors for individual stellar velocities. 

To check consistencies in velocity between 2007 observations and those of 2006, the 2007 
fields were cross correlated with the four 2006 standards as well as the two 2007 standards, 
and the resulting mean stellar velocities and dispersions were compared. Again, due to the 
differing wavelength coverage between 2006 and 2007, the 2006 standards correlated with 
the 2007 targets yielded velocities with an average error higher than the error reported by 
the 2007 standards. We compare the weighted average of the two 2007 velocities with the 
velocities obtained from correlations of the 2006 standards and find, for the 5 fields observed 
in 2007, the individual stellar velocities were offset, on average, by 1 km s~^ with an avergae 
RMS of 2 km s~^, as can be seen in figure [71 Additionally, we re-observed 3 stars in 2007 that 
had been obtained in the 2006 dataset. Two of these 3 stars were found to have velocities in 

2007 that agree with the velocities obtained in 2006, within 4 km s~^, which is less than our 
adopted individual velocity error of 5 km s~^. The third star was found to have a velocity 
difference of 8 km s~^, still in reasonable agreement with the individual errors added in 
quadrature. 

For purposes of model comparison, we follow the lead of Beaulieu et al. (2000) and 
correct all heliocentric velocities obtained for the solar reflex motion using the circular motion 
of the LSR at the Sun as 220 km s~\ and the Sun's peculiar velocity relative to the LSR 
as 16.5 km s"^ toward (/,&) = (53°,25°) (Kerr & Lynden-Bell 1986, Mihalas & Binney 1981, 
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Beaulieu et al. 2000) by 

Vgc = Vhc + 220sin{l)cos{b) + I6.5[sin{b)sin{25) + cos{b)cos{25)cos{l - 53)] (1) 

where Vhc is the reported hehocentric velocity from 'fxcor'. We did not apply this trans- 
formation in Rich et al. (2007a), nor did we apply any additional color cuts (to be discussed 
in section 5), reporting only heliocentric velocities in that study. 

5. Color Selection and Reddening 

In order to assure that our observed stars lie on the RGB, all the field CMDs were 
dereddened using reddening maps provided by Schlegel et al. (1998). For each field, the 
2MASS stars were dereddened as well as our observed targets (figure [8]). The de-reddened 
fields appear reasonable and show a clearly defined "red-edge" of the RGB. Figures [9] and 
[To] show contours of the combined de-reddened 2MASS fields along the b=—4° major axis 
strip and the — 0.4°</<0.0° minor axis strip, along with a greyscale of our de-reddened 
observed stars and isochrones (Marigo et al. 2008). As can be seen, the de- reddened RGB is 
clearly defined, and shows the majority of our sample stars lay within a region of reasonable 
metallicity as expected for bulge stars (e.g. Fulbright, McWilliam, & Rich 1996). In order 
to assure that our observed targets are indeed red giants, we make a color cut omit from our 
sample stars bluer than the [Fe/H]=— 2.0 isochrone (extended upwards), and stars brighter 
than Kr^7 in our calculation of the individual field statistics. Although this reduces our 
observed sample by approximately 4% (see table [211]), to roughly 99 stars per field, this 
selection can be comfortably considered to include members of the RGB of the Galactic 
bulge population. 

For each field, mean velocities and dispersions are calculated directly from a a-clipping 
algorithm, with a 6a cut. None of our observed fields contained any stars with velocities 
deviating by more than 6a, resulting in no stars cut from our color selected sample when 
calculating the statistics. From this sigma-clipping algorithm we obtain for each field a mean 
velocity, velocity dispersion, and calculate the errors in those values as a/y/N and a/\/2N, 
respectively. Of course, the intrinsic dispersion of each field is given by 



N 



intrinsic observed 



Y,[error\v,)]/[2{N-l)] (2) 



1=0 



where the first term is the observed velocity dispersion of a given field and the second term 
represents the uncertainties in the individual stellar radial velocity (v^) measurements, with 
N being the number of stars in that field. For all fields observed, the second term is negligible 
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and thus our observed velocity dispersions can be considered intrinsic to the field. As a 
check of our velocity results, the BRAVA dispersion and mean velocity measurements for 
the Baade's window field (/,&)~(l°,-4°) are (t=112 ± 10 km s~^ and Vmean=-5 ± 14 km s"^ 
and are in good agreement with that of SWC who found values of o"=113 ± 6 km s~^ and 
Vmean=— 4 ± 8 km s~^, which confirms the validity of our measurements. A complete listing 
of our currently observed fields is shown in Table HJ including the data presented in Rich et 
al. (2007) updated with color cuts and galactocentric velocities. 

In order to characterize the kinematic properties of the disk, in comparison to our data, 
we observed a field at (/,&) = (— 30°,— 4°) in 2007. Again, we de-reddened this field, and 
applied the same velocity correction for the solar refiex motion as well as the same color cuts 
applied to our bulge fields to obtain a mean Galactocentric velocity and dispersion (Figure 
[TTi) . Both the velocity dispersion (a = 51 ± 4 km s~^) and mean velocity (—162 km sec^^) 
contrast strongly with our bulge dispersions, which typically exceed 80 km sec~^. Further, 
we note the test applied by Sharpies et al. (1990), in which those stars in Baade's Window 
with I < 11.8 show cr = 71 km s~^, significantly lower than the 113 km s~^ of the full, 
fainter sample. Our disk field, and the study of Sharpies et al. (1990) do not, by themselves, 
prove our sample free of disk contamination. However, the disk field departs so strongly 
from our rotation and dispersion profiles (Figure 12) that we suspect disk contamination is 
not an important issue. Composition measurements and modeling of the CMD, as well as 
comparison of our data to the Besancon galaxy model, will constrain further the level of disk 
contamination. 



6. Results and Analysis 

Figure [12] shows the BRAVA velocities and dispersions compared to the predictions of 
the Zhao (1996) model, in Galactocentric velocity rest frame. The Zhao model is a 3D 
steady-state model which uses a generalized Schwarzchild technique, consisting of orbital 
building blocks within a rapidly rotating bar potential, with corotation at 3.3 kpc, and a 
Miyamoto-Nagai (MN) disk potential (see Zhao 1996 for details). In this model, roughly 
half of the mass of the bar consists of stars on direct regular orbits which form the rotational 
support of the bar, and the remaining half in irregular and retrograde orbits. The data agree 
with the Zhao model dispersion, but shows a significant deviation from the predicted rotation 
curve at large Galactic longitude. Solid body rotation claimed by numerous previous studies 
(Menzies 1990; Izumiura et al. 1995) is not apparent in the 6=— 4° major axis strip; after 
reaching an amplitude of 75 km s~^, the rotation curve fiattens beyond |/|~4°, coresponding 
to a projected distance of ~550 pc from the central minor axis for Ro=8kpc. Furthermore 
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we find that the upper hmit for rotation of the bulge, in Galactocentric rest frame velocity, 
is V/jo(.<100 km s~^ kpc~^ for | / |^4°, and flattens beyond that at a value of roughly 
75 km s~^ (Figure [T^ . Our upper limit for rotation value is significantly larger than the 
previously reported value of VRoi.=70 km s~^ kpc~^ obtained from PNe (Beaulieu et al. 
2000). Interestingly, this velocity flattening is also observed at the same projected distance 
in the proper motion studies of Clarkson et al. (2008, private communication), perhaps 
indicating a change in the density profile. Also plotted in Fig. [12] is the mean velocity and 
dispersion for our observed disk field at (/,&) = (— 30°,— 4°). As can be seen, our disk field 
shows lower dispersion and significantly greater rotation speed, than any of our bulge fields; 
this gives us added reassurance that our sample is not contaminated by foreground/disk 
stars. 

It is also useful to compare BRAVA data to the PNe from the survey of Beaulieu 
et al. (2000), and to follow their example by comparing BRAVA data to other N-body 
models. Two such models are plotted in Figure [121 iii heliocentric velocity. The models 
of Fux (1996) and Sellwood (1993) are both N-body bars formed from initially unstable 
disks, investigating the initial formation of the bulge (see Beaulieu et al. 2000 for details). 
Figure [12] shows that none of the models are satisfactory in predicting both dispersion and 
rotation, suggesting that the data are challenging to fit with either disk-instability formed 
bars or Schwarzchild models with fixed potentials (Zhao 1996). However, there is striking 
agreement in the qualitative inflections of the Fux (1996) dispersion profile and our bulge 
sample, differing only in a relative shift in the scale of the dispersion, while the model 
rotation curve is in excellent agreement with the data. Similarly, the Sellwood (1993) model 
shows qualitative consistency with the observed data; although differing in amplitude, the 
model shows the same flattening seen in the observed rotation curve. Due to the limited 
number of PNe at 6=— 4° and /=|10|°, we accept those in the range of — 8°< b <— 3°. Within 
that longitude/latitude range there remain only 133 PNe for comparison with our data; the 
PNe are then binned so that ~27 PNe reside in each bin (Fig. [T^j). Considering the less 
secure distances and the assignment of population of PNe, the agreement is good and is not 
consistent with solid body rotation for the bulge. 

BRAVA also obtained a minor axis strip of the bulge and has found no evidence of 
minor axis rotation (Figure UM . although more data are needed at positive latitudes (we 
note that the much greater extinction at positive latitudes will make this effort challenging). 
It is also noteworthy that our rotation curve now agrees with the two points observed by 
Minniti et al. 1992 (see discussion in Beaulieau et al. 2000), but we emphasize that for 
{l,b) = {12° , 3°) Minniti et al. (1995) describe the field as "heavily contaminated by disc 
stars" and we consider any agreement between this field and BRAVA to be fortuitous; the 
connection between the bulge and fields with \l\ > 10° demands more investigation (see 
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Figure 7 of Beaulieau et al. 2000). Indeed, inspection of figure 7 of Beaulieu shows a 
break in velocity and dispersion beyond |/| > 10°, suggesting two different populations 
being sampled. Reassuringly, our disk field at (/,&) = (— 30°,— 4°) matches the PNe data quite 
well and illustrates the importance of gathering more data to investigate the kinematics of 
fields between 10°< |1| <30°. In our opinion, our 2MASS selection method is better able 
than optical colors to defeat the reddening and to separate cleanly the bulge giants from 
foreground red giant and clump stars. At (/,6) = (12°,3°), it would be critical to have sampled 
adjacent fields and to have established a connection between this field and the higher surface 
brightness, inner bulge fields. At this low latitude, disk contamination becomes a concern in 
any optically selected sample. 

To check for metallicity and/or magnitude bias in our survey, we have divided our 
de-reddened sample by magnitude (K<8.7 and K>8.7) as well as color (see Figure [T^ . 
Dividing our sample by magnitude seems to have little effect upon the kinematics of the 
major axis strip, further reinforcing our belief that our sample does not experience significant 
contamination from other foreground stellar populations, to first order. Dividing our sample 
by color shows little effect as well, however, there appears to be a general trend indicating 
that the more metal rich giants show a lower dispersion in each individual field. Of course, 
by dividing our sample, we reduce the statistical significance of this apparent trend. We 
plan to measure the spectroscopic metallicities of our sample. 

In 2006, hints of kinematic substructure in the velocity distributions of several fields 
were noted and reported in Rich et al. (2007b). We now present the Galactocentric velocity 
distributions of all of our current bulge fields in figures [16] and [T71 As can be seen, ~ 22% 
of the observed bulge fields show 2.5(T deviations, as estimated by eye, from a Gaussian 
constructed from a calculation of mean velocity and dispersion in each field. Indeed, the 
choice of bin width determines the amount of "structure" seen and at a bin width of 25 
km s~^, we are perhaps overbinning our data. The choice of bin width is driven by the ~ 10 
km s~^ dispersion expected for cold components like a dissolving cluster or dwarf galaxy. 
Therefore, we believe it is still useful to display the data in this fashion in order to easily pick 
out features that may justify follow-up observations. Although the third (skew) and fourth 
(kurtosis) moments of the assumed Gaussian distributions calculated for each individual field 
was not significant, those fields with apparent deviation from a normal distribution warrant 
a closer inspection. Possible substructure might arise from disrupted satellites or stars in 
unique orbit families but the relatively small number of stars observed in each field require 
subsequent observations to increase our sample size. Our 2007 observations were aimed at 
investigating two of these fields to determine if the 2.5cr deviations from the Gaussian in 
select 2006 fields (see Figure [18]) were real or merely statistical in nature. Two approaches 
were used to test these possible kinematic clumps. For the field (/,6) = (— 8°,— 4°), we shifted 
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by half a degree in b above and below the original field, in order to detect a possible spike at 
the same velocity seen in the initial observation. For the field (/,6) = (4°,— 4°), we obtained two 
separate exposures at the same (l,b), but selecting different stellar targets from the 2MASS 
catalog. As can be seen in Figure [TSl the initial spikes turned out to be statistically not 
significant. However, we intend to re-observe other fields which exhibit similar structures. 
Simulations by Reitzel et al. (2007) show that such deviations from a normal distribution 
are not unexpected in random draws, and preliminary examination suggested that roughly 
16% of the apparent substructures from our data set might be real. A full statistical analysis 
between models and our data will be able to determine if we are seeing more "spikes" than 
are expected from random noise, but we will be unable to determine which one is real without 
gathering more data for each field. It is important to note that a dissolving system may have 
a range in abundance, so that the addition of abundance would not necessarily increase the 
significance of a detection. 

We now turn to the coadded total velocity distribution. We have constructed the minor 
axis Galactocentric rest frame velocity distribution, giving N mean = —2 ± 4 km s"'^ and a = 
117 ±3 km s~^ for ~820 stars, as well as having coadded the entire sample in Galactocentric 
coordinates, yielding Vmean = 1 ± 2 km s~^ and a = 116 ± 1 km s~^ for ~3200 stars (Figure 
[T9l) . In both cases, the distribution appears Gaussian, with no apparent deviation from a 
normal distribution greater than f»1.5(T, as estimated by eye. For the entire coadded bulge 
sample, when fit by a x^ minimization technique, we find a reduced x^ = 1.29, with a mean 
of 0.6 ± 2.4 km s~^ and a = 122 ± 2 km s~^, consistent with our calculations of the first and 
second moments of the distribution. Furthermore, our distribution has a (negligible) skew 
of 0.05 ± 0.04 km s~^ and a kurtosis of —0.36 ± 0.09 km s~^, suggesting a slightly platykurtic 
(fiattened) distribution. The negligible skew is in contrast to the measured distributions 
of SWC (1990; Figure 9) and Minniti (1992; Figure 1), which show strongly fiattened or 
significantly skewed characteristics, indicating perhaps contaminated samples. The lack of 
such characteristics in our distributions bolsters the argument that our samples are not 
contaminated by cold components (disk) or hot components (halo). 

Furthermore, we have investigated characterizing our co-added samples using the Gauss- 
Hermite series as defined by van der Marel & Franx (1993), which yield mean velocity, dis- 
persion, as well as two other dimensionless parameters that measure the asymmetric (h_3, 
analogous to skew) and symmetric (h_4, analogous to kurtosis) deviations of a distribution 
from Gaussian (Figure [T9|) . According to van der Marel & Franx (1993), this method mini- 
mizes correlations between the parameters of the fit, and allows for a more sensitive measure 
of the deviation from a Gaussian line profile. A least squares fitting to each summed distri- 
bution yields V^ean = -2.6 ± 5.4 km s-\ a = 122.2 ± 5.3 km s-\ h_3= 0.038 ± 0.031, and 
h_4= 0.002 ±0.031 for the minor axis, and Vniean= 0.5 ±2.6 km s-\ a = 121.9 ±2.5 km s~\ 
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h_3= 0.011 ± 0.015, and h_4= —0.014 ± 0.015 for the entire bulge sample. These values are 
in excellent agreement with our more standard calculation of the moments of an assumed 
Gaussian distribution for the Galactic bulge. Figure EO shows our longitude- velocity plot 
for the major axis strip. As can be seen, there is no evidence of a cold, disk component 
in our sample which would manifest itself as a linear trend. There is no indication of a 
departure from a Gaussian single population in any field. The overall trapezoidal appear- 
ance of the outer envelope of the distribution is reminiscent of the gas dynamics observed 
by Liszt & Burton (1980) within 2 kpc; the gas motions were modeled as arising in a tilted 
bar potential. Our repeat tests with gaussian fitting, examination of the first four moments 
of the distribution, as well as the fits using the Gauss- Hermite series are all consistent with 
a purely Gaussian distribution, with no indication of contamination from any foreground 
stellar population. Furthermore, it is noteworthy that we discovered no stars with velocities 
greater than ±4o" in any of our fields. 



6.1. The Bulge on the (V^a^/a) vs. e Plot 

Armed with the velocity dispersion and a rotation speed, it is now possible to use the 
BRAVA measurements to place the Galactic bulge on the Kormendy & Kennicutt (2004) 
plot of (Vmax/o') VS. € (Figure!^. One interprets the ratio (Vmax/c") as a ratio of ordered 
to random kinetic energy and e as the apparent flattening of the bulge. Placement of galax- 
ies on this plot represent a method of interpreting the formation of bulge systems; bulges 
formed via secular proceses (pseudo-bulges) lie on or above the oblate line and are more ro- 
tation dominated than classical bulges and ellipticals which are formed via mergers. A very 
thorough discussion can be found in the review of Kormendy & Kennicutt (2004), and shall 
not be repeated here. We use the Weiland et al. (1994) minor-to-major axis ratio of ~ 0.6 
to estimate e = 1— axial ratio= 0.4. We adopt o"=116 km s~^ for the velocity dispersion, 
and 75 km s~^ for the rotation speed, giving (Vmax/o")=0.64±0.5. This places our Milky 
Way near the well-known edge-on spiral NGC 4565 with its peanut-shaped bulge, and near 
the oblate line, but below the most rapidly rotating bars and pseudobulges. Kormendy & 
Illingworth (1982) note that NGC 4565 not only has a peanut-shaped bulge, but exhibits 
cylindrical rotation (no decrease in rotation speed as a function of vertical distance above 
the plane), and therefore (Vmax/o') underestimates the dynamical importance of rotational 
kinetic energy compared to the ellipsoidal bulges. While residing among the classical bulges 
in this plot, NGC 4565 and the Milky Way share the classic peanut bulge morphology while 
still lying well below the oblate line. Kormendy & Kennicutt (2004) note that a peanut- 
shaped bulge is sufficient to classify a galaxy has having a pseudobulge, so it may be no 
surprise that the Galactic bulge would lie beneath the oblate line should it indeed exhibit 
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cylindrical rotation similar to NGC 4565. While the stellar population of the bulge is old 
(see e.g. Ortolani et al. 1995) a pseudobulge can be old in the Kormendy & Kennicutt 
(2004) definition. The measurement of cylindrical rotation in the M giant population would 
strengthen substantially the analogy between the Milky Way's bulge and that of NGC 4565; 
indeed, we have observed this in the —8° latitude fields (Howard et al. 2008, in prep). We 
defer consideration of the Milky Way's putative pseudobulge for this future work. 



7. Conclusions 

We present our methods of selecting M giant probes for the Bulge Radial Velocity Assay. 
Our sample as reported in this paper has nearly 3300 stars, and is defined clearly to cover 
the full range of abundance on the first red giant branch. We report our rotation curve and 
velocity dispersion profiles, which agree more closely with the self consistent model of Zhao 
(1996) compared to N-body models in both rotation curve and dispersion. 

The reliability of our year-to-year measurements, along with our velocity precision, allow 
us to search for evidence of cold streams in the bulge resulting from recent merger events, 
analogous to the Sagittarius Dwarf spheroidal stream. We find no evidence for hidden 
dynamically cold, or hot, subcomponents, either in the line of sight velocity distribution of 
the minor axis fields, or in the 1-v plot. Two candidate cold stream features are investigated 
in follow-up studies that increase the sample size by a factor of three. The features were not 
confirmed, and illustrate the importance of sample size when searching for such features. 

We find a departure from solid body rotation at a projected distance of ?ii0.6 kpc for 
Ro=8 kpc, in agreement with the proper motion study of Clarkson et al. (private communi- 
cation). The break from solid body rotation occurs as the density of the bulge is dropping 
and may signify a transition from bar to keplerian dynamics. The bulge appears to extend 
to / = ±10°, beyond which a colder, more rapidly rotating, disk-like population appears to 
dominate. However, it is noted that there is a lack of data between 10°<| / |<30°, which 
prevents one from determining the relationship between our bulge sample, and the outer disk 
fields. Our selection technique has also allowed us to determine that the the bulge population 
is normally distributed, with no significant deviation from Gaussian when considering the 
entire co-added sample of ~ 3200 stars. 

We have divided our sample by magnitude and color, and find no significant differences 
in dynamics. We believe that our selection technique is robust and is yielding a bulge/bar 
dominated sample. The rotation curve permits us to place the Milky Way's bulge on the 
(Vmax/o') VS. e plot; it is close to the peanut-shaped bulge of NGC 4565 in this diagram. 
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Further observations at the b= —8° latitude fields will allow us to determine if the Milky 
Way bulge rotates cylindrically, an indication of whether the Galactic bulge is a pseudobulge, 
or a more classically evolved bulge. 
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Fig. 1. — A simple schematic of the bulge/bar orientation with respect to the Sun/Galactic 
center hne of sight. Exact values of the angle are uncertain but range from ~15° to ~35° 
(Alcock et al. 2000; Gerhard et al. 2002). On this diagram, positive Galactic latitude {b) 
points up out of the page. 
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Fig. 2. — Observed BRAVA fields, up to April 2007, overplotted on the COBE 2/im image 
(Launhardt et al. 2002). Grey circles represent fields observed in 2005, black circles are 
fields observed in 2006, and white circles are fields observed in early 2007. The size of the 
circles corresponds to the 40" field of view of the instrument. Dashed circles represent fields 
observed over multiple years. 1° corresponds to 140 pc for a distance of 8 kpc; the h= —4° 
strip lies roughly 550 pc south of the Galactic center. The Galactic plane is avoided due to 
heavy extinction. 
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Fig. 3.— Color-Magnitude Diagram of 2MASS targets in the field at (l,b) = (4°,-4°). Filled 
circles indicate stars for which we have spectra. Included are reddened isochrones (Girardi 
et al. 2002) for [Fe/H]= —1.3 and —0.5. The reddening vector corresponds to E(J-K)=0.33 
from the Schlegel et al. (1998) map. The parallelogram indicates our selection region; the 
blue cutoff rejects many objects that are closer than the bulge, which have lower reddening 
and are brighter than the red giant branch. 
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Fig. 4. — A sample of the spectra observed from each year of observations, 2005 at top, 
moving down to 2007 at bottom. The lower panel is a blue star observed in 2006 to illustrate 
where there are areas of significant atmospheric contamination. The hash marks below the 
blue star mark regions used for cross correlation. Evident in the stellar spectra are the TiO 
bands at ~7050A , as well as the Telluric A band at ~7600A. In the 2006/2007 spectra the 
Call triplet is seen at ~8498A,~8542A, and, beyond the spectral range covered in 2006 but 
observed in 2007, the third line at ~8663A. 
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Fig. 5. — A typical cross correlation peak as reported by the IRAF cross- correlation routine 
'fxcor'. As can be seen in this example from 2006, the correlation peaks are quite strong, 
and yield Tonry-Davis R (TDR) parameters of 10-30, yielding velocity errors between 5-15 
km s~^, depending on the instrument set-up used. 
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Fig. 6. — Stellar radial velocity offset between 2005 and 2006 observations for field (l,b) = (6°,- 
4°). Our 2005/2006 individual stellar velocities show agreement to within 5kni s~^. We 
therefore adopt 5 km s~^ as our individdual stellar velocity error. Since the dispersion we 
see in our fields is over an order of magnitude larger than out individual stellar velocities, we 
consider these velocity errors to be negligable. The error bar in the top left of the figure shows 
the horizontal 2006 errors as reported by 'fxcor' as well as the 2005/2006 errors reported by 
fxcor, added in quadrature. 
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Fig. 7. — Comparison of velocities for the 2007 observations. For the 5 fields observed in 
2007A, velocities were obtained from the four 2006 standards and the two 2007 standards. 
The vertical axis are the difference in velocities and each panel shows the average offset and 
RMS scatter of the offset. As can be seen, our velocities appear to be, on average, consistent 
to within our adopted error of 5 km s~^. 
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Fig. 8. — Sample grayscale CMDs of four of the 35 BRAVA survey fields. The upper panels 
of each figure show all 2MASS targets obtained in each field, with BRAVA targets marked 
with points. The lower panels show the same ^M^dS'^' stars and BRAVA targets, individually 
de-reddened using Schelgel et al. (1998) reddening maps. As can be seen, some fields show 
observed stars that appear too blue in color to be safely assured of their bulge membership. 
These blue stars are omitted when calculating the field statistics (see figures 1^ . [TUl) . 
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Fig. 9. — Left: Color-Magnitude contours of the 2MASS catalog (~440,000 stars) for our 
fields along 6=— 4° major axis strip. Overlayed in greyscale are our observed BRAVA targets, 
representing 2612 stars, and isochrones for a 12 Gyr population at a distance modulus of 
14.47 mag. The isochrones, starting on left, are for [Fe/H]= —2.0,— 1.3,— 0.5, and +0.2 
(Marigo et al. 2008). Right: Same as left, but showing the BRAVA targets remaining after 
trimming stars bluer than the [Fe/H]= —2.0 isochrone and brighter than K~7.4. Out of the 
original 2612 targets observed along this strip, 2505 survive the cut. These 2505 targets are 
used in calculation of individual field statistics. 



-30- 





Fig. 10.— Left: Contours of the 2MASS catalog (~290,000 stars) for our fields along 
— 0.4°</<0.0° minor axis strip. Overlayed in greyscale are our observed BRAVA targets, 
representing ~850 stars, and isochrones for a 12 Gyr population at a distance modulus of 
14.47 mag. The strong effects of reddening and the uncertaintly of the Schelgel maps at 
h< |5|° can be seen in the broadening of the RGB, as compared to fig. [9] The isochrones, 
starting on left, are for [Fe/H]= -2.0,-1.3,-0.5, and +0.2 (Marigo et al. 2008). Right: 
Same as left, but showing the targets after trimming stars bluer than the [Fe/H]= —2.0 
isochrone and brighter than K 7.4. Out of the original 850 BRAVA targets observed along 
this strip, 832 survive the cut. These 832 targets are used in calculation of individual field 
statistics. 
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Fig. 11.— Left: Greyscale CMD of our disk field at (/,&) = (-30°,-4°), with tfie observed 
targets marked as points. Tlie same color cuts and dereddening applied to our bulge sample 
set are applied here. The isochrones are the same as those described in fig. [9l Right: 
Kinematic data of our disk field in the Galactocentric rest frame. Binsize is 25 km s^^. 
Comparison with Table H] and Figure [12] show this disk field to have a lower dipsersion and 
mean velocity that the fields in the bulge. 
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Fig. 12. — Left: Velocity dispersion profile and rotation curve of bulge major axis, b=—4°, 
in galactocentric velocity; The solid line indicates the model of Zhao (1996). Zhao's model 
satisfactorilly predicts the dispersion, but shows a more rapid rotation that what is observed. 
Our disk field at (/,&) = (— 30°,— 4°) has been plotted to show the lower dispersion from a 
disk field, as expected. Right: Major axis velocity dispersion profile and rotation curve, 
in Heliocentric velocity. Also indicated are the PNe data of Beaulieau et al. (2000) (" + " 
symbols) and the models of Fux (1997) (dotted fine), and Sellwood (1993) (dashed line). Our 
velocity dispersion are higher than those of the PNe sample, while our rotation is consistent 
with the PNe. As a homogenous population, the PNe may suffer some disk or forgeground 
contamination. For the first time, our densely sampled fields follow the qualitative infiections 
of the N-body dispersion and rotation profiles of Fux (1996) and Sellwood (1993) quite well. 
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Fig. 13. — Best fit line of tlie central 8° of the Galactic bulge. The line shows an upper limit 
of the slope of VRot. = 100 km s~^ kpc^^. The break from solid body rotation is evident at 

I /I -4°. 
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Fig. 14. — Comparison of data with the Zhao (1996) model of the Galactic bulge minor axis 
(6=0°,— 0.4°< / <0°). Top- The mean velocity is consistent with the model in that it shows 
no appreciable minor axis rotation. Bottom-The velocity dispersion is consistent with the 
trend of the model at negative latitudes, although more data is needed at positive latitudes. 
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Fig. 15. — Kinematics of the bulge major axis, separated by color /magnitude. The top 
plot shows the total de-reddened 2MASS CMD of the bulge in countour, with our observed 
targets in greyscale. The dashed lines represent our cuts for metallicity and brightness. As 
can be seen, there appears to be no bias in our kinematics based on brightness. The color 
cut offers a hint that the red (metal rich) population has a slightly lower velocity dispersion 
in some fields; followup with larger samples and actual metallicity measurement is required 
to confirm or refute this. 
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Fig. 16. — Presentation of all bulge field Galactocentric velocity distributions. Overlaid on 
each plot is a Gaussian with parameters derived from the IDL cr-clipping algorithm 'meanclip' 
used to calculate field statistics. As can be seen, several fields display a departure from a 
normal distribution, most notably fields 2,3,7,8, and 12. However, we stress that these 
departures can be a result of choice of bin size. Fields 3 and 12 were chosen for follow-up 
observation in 2007; the deviation in those fields were not statistically significant (see fig. 
[18] for details). Not shown are the velocity distributions of our disk field (fig. [TTl) . or field 
P4-2005 which has the same stars observed as field P4-2006. Bin size is 25 km s~^. 
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Fig. 17. — Bulge line of sight velocity distributions cont'd, (see fig. [16 
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Fig. 18. — Two candidate stream fields from 2006 were re-observed in 2007. For each field 
location, the individual fields and statistics are shown in the first three panels, and the final 
summed field is shown in the final panel. As can be seen, our initial detections of "spikes" 
turned out to be not significant, and emphasizes the importance of sample size when looking 
for cold stream features. Bin size is 25 km s""*^. 
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Fig. 19. — Left: All minor axis fields summed, consisting of 822 stars. Right: All bulge 
fields (both minor and major axes) summed, consisting of 3202 stars. Both figures are in 
Galactocentric coordinates with a Gaussian with Vmean and a derived from the IDL from 
a (j-clipping algorithm, 'meanclip', overlayed (solid line). Also plotted is the fit of the 
Gauss- Hermite series (dotted line) as described in Section 6. As can be seen, both curves 
are consistent with each other, and suggest our bulge sample consists of a homogeneous, 
normally distributed, stellar population. Furthermore, measurements of the 3rd and 4th 
moments of the Gaussian distribution as well and the h_3 and h_4 parameters of the Gauss- 
Hermite series suggest no significant symmetric or assymetric deviations from a Gaussian 
distribution. 
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Fig. 20. — Longitude- Velocity (1-v) plot for the entire bulge sample, smoothed to 1° in 
longitude and 10 km s~^ in galactocentric velocity. The higher density of stars at 1=0° is 
due to the inclusion of the minor axis fields, while the higher density at 1=— 8° and 1=4° are 
due to our follow-up observations of those two fields in 2007 (see fig. [TS]). As can be seen, 
there is no evidence of a cold, disk component in our sample which would manifest itself as 
a linear trend. We also do not see any evidence of a hot component that would indicate 
presence of a halo/spheroid population. The roughly trapeoidal envelope of this distribution 
is reminiscent of that seen for the gas dynamics at R<2 kpc (Liszt & Burton 1980). 
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Fig. 21. — (Vmaaj/o") plot from Kormendy & Kennicutt (2004) with the Galactic bulge indi- 
cated (black cross). The MW bulge lies under the oblate supported line and less rotationally 
supported than the pseudobulges (filled symbols) but is similar to classical bulges (open 
symbols). NGC 4565 is an edge-on spiral, similar to the Milky Way, in that is exhibits a 
box-shaped bulge. 
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Table 1. BRAVA Observation Details 



BRAVA field 


Date of Obs. 


Pointing 


Center 




Exposure 


Grating 


Central A 


Slit Mask 




(UT) 


(J2000) 






Angle 


(A) 




fa4 


Jul 30 2005 


18'' 00"" 25.02= 


-31°58'53' 


.1 


Ix300s,lx900s 


18.206 


~7600A 


_ 


fe4 


Jul 30 2005 


17''51'"43.92= 


-35°19'00' 


.1 


3x300s 


18.206 


~7600A 


- 


fm4 


Jul 30 2005 


18''08'" 17.54= 


-28°21'29' 


.5 


3x300s 


18.206 


~7600A 


- 


fo4 


Jul 31 2005 


18''12'"31.04= 


-26°36'49' 


.4 


2x420s 


18.206 


~7600A 


- 


fp4-2005t 


Jul 31 2005 


18''14'"40.35= 


-25°42'48' 


. 1 


3120s* 


18.206 


~7600A 


- 


fMA 


Aug 


01 2005 


17'' 53^27.22= 


-29°58'22' 


.4 


2x900s,lx600s 


18.206 


~7600A 


- 


fMk 


Aug 01 2005 


17'' 28^31.48= 


-26°31'40' 


.2 


2x900s 


18.206 


~7600A 


- 


fl 


May 


16 2006 


17''36'"27.20= 


-39°30'23' 


.8 


3x900s 


18.732 


~7800A 


200Atm 


fl6 


May 


16 2006 


18'' 22^52. 16= 


-22°10'36' 


.1 


3x900s 


18.732 


~7800A 


200Atm 


f8 


May 


16 2006 


17'' 56'" 53.48= 


-32°41'07' 


.7 


3x900s 


18.732 


~7800A 


200^m 


fl2 


May 


17 2006 


18'' lO"* 19.68= 


-27°28'48' 


.1 


2x600s,lx900s 


18.732 


~7800A 


200^m 


fl5 


May 


17 2006 


18''20'"47.74= 


-23°03'57' 


.9 


2x300s,lx600s 


18.732 


~7800A 


200Atm 


f3 


May 


17 2006 


17''41"' 55.61= 


-37°50'19' 


.3 


3x600s,lx300s 


18.732 


~7800A 


200^m 


f5 


May 


17 2006 


17''47'"03.31= 


-36°08'10' 


.3 


3x600s 


18.732 


~7800A 


200Atm 


fswc 


May 


17 2006 


18''03"' 44.62= 


-30°01'50' 


.0 


3x600s 


18.732 


~7800A 


200/im 


flOl 


May 


18 2006 


17''49"'43.50= 


-29°56'20' 


.9 


2x600s,2xl200s 


18.732 


~7800A 


200Atm 


fl3 


May 


18 2006 


18''16"'39.90= 


-24°50'11' 


.1 


3x600s 


18.732 


~7800A 


200^m 


fl4 


May 


18 2006 


18''18"'46.62= 


-23°57'49' 


.3 


3x600s 


18.732 


~7800A 


200Atm 


f4 


May 


18 2006 


17'' 44^32.07= 


-36°58'46' 


.1 


3x600s 


18.732 


~7800A 


200Atm 


f7 


May 


18 2006 


17'' 54"' 26. 79= 


-33°33'12' 


.1 


3x600s 


18.732 


~7800A 


200^m 


fll 


May 


19 2006 


18''06"'03.16= 


-29°13'30' 


.6 


3x600s 


18.732 


~7800A 


200Atm 


f2 


May 


19 2006 


17''39"'12.63= 


-38°41'39' 


^7 
. 1 


3x600s 


18.732 


~7800A 


200Atm 


f6 


May 


19 2006 


17''52"'04.82= 


-34°25'09' 


.3 


3x600s 


18.732 


~7800A 


200Atm 


f9 


May 


19 2006 


18''01"'30.41= 


-30°58'26' 


.1 


3x600s 


18.732 


~7800A 


200^m 


fp4-2006t 


May 


19 2006 


18'' 14^40.35= 


-25°42'48' 


. 1 


3x600s 


18.732 


~7800A 


200Atm 


f27 


May 


20 2006 


18''05"'38.05= 


-31°26'59' 


.1 


3x600s 


18.732 


~7800A 


200^m 


f63 


May 


20 2006 


18''08"'42.54= 


-31°51'07' 


.0 


3x600s 


18.732 


~7800A 


200^m 


f66 


May 


20 2006 


17''56"'46.67= 


-30°43'32' 


.1 


4x600s 


18.732 


~7800A 


200Atm 


f90 


May 


20 2006 


17''59"'33.84= 


-30°43'21' 


.1 


3x600s 


18.732 


~7800A 


200Atm 


fSOSO 


Apr 07 2007 


16''23"'00.60= 


-55°17'29' 


.8 


3x600s 


18.985 


~7900A 


200^m 


fa 


Apr 25 2007 


17''41"' 14.27= 


-38°00'06' 


.6 


5x600s 


18.983 


~7900A 


200Atm 


fb 


Apr 25 2007 


18''08"'25.38= 


-27°13'55' 


•y 
. 1 


4x600s 


18.983 


~7900A 


200Mm 


fc 


Apr 


25 2007 


17''41"'35.59= 


-37°56'12' 


.6 


3x600s 


18.983 


~7900A 


200^m 


fd 


Apr 


25 2007 


18''12"'29.34= 


-27°42'40' 


.6 


3x600s 


18.983 


~7900A 


200Mm 



tpield P4 observed with the same configuration in 2005/2006 
*3120s=lx300s,lx420s,lx600s,2x900s 
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Table 2. Number of Reliable Velocity Measurements per Field 



BRAVA field 



Pointing Center 

(l°,b°) 



Target Fibers in 
Configuration 



Sky Fibers in 
Configuration 



Reliable Velocities^ 



fa4 


-1.01, -4.29 


110 


20 


fe4 


-4.81, -4.40 


105 


25 


fm4 


2.99, -4.01 


111 


21 


fo4 


4.98, -4.00 


109 


21 


fp4-2005 


6.00, -4.00 


106 


20 


fMA 


-0.01, -2.00 


111 


20 


fMk 


0.01, 4.48 


100 


26 


fl 


-9.98, -3.98 


96 


26 


fl6 


10.02, -3.99 


107 


21 


f8 


-1.99, -3.99 


109 


19 


fl2 


3.98, -3.99 


107 


24 


fl5 


9.00, -3.99 


105 


22 


f3 


-8.00, -4.00 


103 


25 


f5 


-6.00, -3.99 


106 


22 


fswc 


1.05, -3.96 


74 


20 


flOl 


-0.40, -1.28 


100 


31 


fl3 


6.99, -3.98 


105 


22 


fl4 


7.99, -4.00 


108 


22 


f4 


-6.99, -4.00 


104 


25 


n 


-3.00, -3.98 


112 


19 


fll 


2.00, -4.00 


110 


19 


f2 


-9.01, -4.00 


102 


24 


f6 


-4.00, -4.00 


109 


22 


f9 


-0.01, -4.00 


110 


20 


fp4-2006 


6.00, -4.00 


109 


21 


£27 


0.00, -5.00 


106 


23 


f63 


-0.04, -5.77 


106 


22 


f66 


-0.30, -2.99 


109 


20 


f90 


0.00, -3.51 


112 


20 


fSOSO 


-29.99, -3.96 


100 


20 


fa 


-8.21, -3.97 


108 


17 


fb 


4.00, -3.50 


112 


15 


fc 


-8.11, -4.00 


105 


21 


fd 


4.01, -4.52 


111 


17 



109 
101 
108 
109 
101 
108 
95 
95 
106 
109 
105 
101 
102 
105 
72 
97 
104 
107 
103 
109 
107 
101 
108 
110 
108 
105 
106 
107 
112 
96 
100 
111 
102 
109 



t Number of reliable velocities prior to color cuts. See section 5 for details 
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Table 3. 2006/2007 Velocity Standards 



Indentifier 


RA 


DEC 


Spectral Type 


Velocity 


Year Observed 




(J2000) 


(J2000) 




(km s-i) 




HD 203638 


2l'^24"09.593^ 


-20°51'06".73 


KOHI 


22.0 


2006 


HD 177017 


19'^03™43.756« 


-22°42'43".00 


M7ni 


42 


2006 


HD 207076 


21^46^31.849^ 


-02°12'45".92 


M8IIIv 


-37.2 


2006 


HD 218541 


23'"09™05.561^ 


-30°08'02".15 


M6ni 


30.0 


2006 


HD 134140 


15'^08™57.514« 


-26°29'50".18 


Mini 


34.8 


2007 


HD 146051 


16'^14™20.739^ 


-03°41'39".56 


M0.5HI 


-19.9 


2007 
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Table 4. BRAVA Rotation and Dispersion Results 



BRAVA field l,b <Vhc > <yGC > err(<V>) a crr(cr) N* (N**) 

(degrees) (km s"-"-) (km s"-'-) (km s"-') (km s"-"-) (km s"-"-) 

6.1 51.3 
9.3 84.1 

8.3 79.7 

9.2 86.8 
8.7 83.2 

8.4 81.1 
5.1 84.2 

8.5 85.7 
8.9 88.6 

9.0 90.3 

9.3 96.5 

10.7 111.9 
11.0 112.8 

11.4 119.1 

15.0 135.4 
11.9 121.6 

10.5 108.1 

10.8 112.8 

12.1 125.8 

9.7 99.0 
11.8 124.6 

11.6 112.9 
14.3 111.5 

11.3 116.3 
10.1 104.7 

10.4 106.1 
10.3 106.4 
9.5 98.8 

5.8 103.5 

9.1 94.8 
8.5 88.2 
9.1 91.6 
9.5 91.7 
8.5 88.2 
8.1 77.8 
8.8 87.0 

t Summed result of previous 3 fields. For plotting, these values are used for this {l,b). See also Fig. 1181 

^Number of velocities after color cut. This value is used to calculate individual field statistics and err(<V>) 
and err((T). 

**Numbcr of reliable velocities before color cut. This value is the same as reported in Table [2] 



f3030 


-30.00, 


, -3.98 


-53.3 


fl 


-9.97, 


-3.99 


-66.9 


f2 


-9.01, 


-4.01 


-35.5 


fa 


-8.20, 


-3.98 


-68.5 


fc 


-8.09, 


-4.00 


-64.0 


f3 


-7.97, 


-4.00 


-45.9 


f3-allt 


-8.09, 


-3.99 


-59.4 


f4 


-6.97, 


-3.99 


-62.6 


f5 


-6.00, 


-3.99 


-44.6 


fe4 


-4.81, 


-4.38 


-25.8 


f6 


-3.99, 


-3.98 


-44.7 


f7 


-3.00, 


-3.99 


-46.7 


fS 


-2.00, 


-3.98 


-23.2 


fa4 


-1.00, 


-4.28 


-9.9 


flOl 


-0.39, 


-1.29 


-1.5 


166 


-0.30, 


-2.99 


-27.8 


f63 


-0.02, 


-5.78 


-11.4 


f9 


-0.02, 


-3.99 


-14.9 


fMA 


-0.01, 


-2.00 


-16.8 


f27 


-0.01, 


-5.00 


-11.9 


f90 


0.00, 


-3.51 


5.4 


fMk 


0.00, 


4.50 


-1.8 


fswc 


1.05, 


-3.96 


-4.7 


fll 


2.00, 


-4.00 


30.5 


fm4 


3.01, 


-4.01 


37.1 


fl2 


3.99, 


-3.99 


29.2 


fb 


4.00, 


-3.50 


28.3 


fd 


4.01, 


-4.50 


26.0 


fl2-allt 


4.00, 


-4.00 


27.8 


fo4 


4.99, 


-3.99 


31.8 


P4-2006 


6.01, 


-3.99 


49.3 


P4-2005 


6.02, 


-4.00 


44.9 


fl3 


6.99, 


-3.98 


35.7 


fl4 


7.99, 


-3.99 


44.0 


fl5 


9.01, 


-3.98 


22.2 


fl6 


10.00, 


-3.99 


35.7 



161.7 


-98.6 


-63.4 


-93.1 


-88.1 


-69.6 


-83.6 


-82.2 


-60.4 


-36.8 


-52.3 


-50.4 


-22.8 


-5.5 


5.8 


-20.4 


-3.2 


-6.5 


-8.1 


-3.6 


13.9 


7.7 


8.0 


47.0 


57.8 


53.7 


53.0 


50.6 


52.4 


60.4 


82.0 


77.7 


72.3 


84.6 


66.8 


84.2 



4.3 


70(96) 


6.6 


81(95) 


5.8 


93(101) 


6.5 


89(100) 


6.2 


91(102) 


5.9 


94(102) 


3.6 


272(300) 


6.0 


102(103) 


6.3 


100(105) 


6.4 


101(101) 


6.6 


108(108) 


7.6 


109(109) 


7.7 


106(109) 


8.1 


109(109) 


10.6 


82(97) 


8.4 


105(107) 


7.4 


106(106) 


7.6 


110(110) 


8.6 


102(108) 


6.8 


105(105) 


8.4 


111(112) 


8.2 


95(95) 


10.1 


61(72) 


8.0 


106(107) 


7.1 


108(108) 


7.4 


104(105) 


7.3 


106(111) 


6.7 


108(109) 


4.1 


318(325) 


6.4 


109(109) 


6.0 


108(108) 


6.4 


101(101) 


6.7 


93(104) 


6.0 


107(107) 


5.7 


92(101) 


6.2 


97(106) 



